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Sakai, Motome, and Imada Reply: In Ref. U], we have clar- 
ified the structure of poles and zeros of Green's function for 
the 2D Hubbard model in the full momentum-energy space by 
using a cluster extension of the dynamical mean-field theory, 
and have elucidated that the interference between poles and 
zeros plays a crucial role in the doped Mott insulator. The 
study provides a comprehensive understanding of a number 
of unconventional features of Mott physics such as hole pock- 
ets, pseudogap, Lifshitz transitions, non-Fermi liquids, Fermi 
arcs, and a spectral weight transfer over the Mott gap. 

In the preceding Comment yl] on our Letter |l[], Phillips 
claims about only one issue confined to the spectral weight 
transfer, among our wide-ranging clarifications. The criticism 
is that we consider an excess spectral weight just above the 
chemical potential in doped states as a non-trivial finding even 
though it has already been discussed in the prior works yl 0, 
5]. Phillips also criticizes that our argument referring to the 
double occupancy rid is not compatible with the perturbation 
theory [6] and offers no underlying physical ground. 

The former comment has nothing to do with our contribu- 
tion from a careful analysis of updated numerical results. As 
clearly stated in our Letter yj], our emphasis is not on the 
excess weight itself but on its sharp rise by tiny doping. As 
shown in Fig. 3(b) of Ref. |Q]], our numerical data indicate 
that the spectral weig ht Wi (A in Ref. 01) increases quite 
rapidly as a function of doping 5 in the vicinity of the Mott 
transition. Phillips also criticized that our result is identical to 
that in Fig. 2 of Ref. |5], but this is totally incorrect: The re- 
sult in Ref. [5] is for ID system, and furthermore, such rapid 
increase near 6 = is not observed even though there is an ad- 
ditional weight to the static component proportional to 8. The 
very rapid, non-linear behavior found in our result is quite un- 
expected, and at least, to our best knowledge, has not been 
anticipated in previous studies. 

As to the latter criticism on the incompatibility of ours with 
the perturbation result with respect to t/U, it is indeed difficult 
to expect that such abrupt change in W\ is described by per- 
turbative arguments. Although Phillips criticized the incom- 
patibility, the first-order perturbation, which is reproduced in 
Eq. (1) in Ref. H, trivially fails as S — » 0, because the ki- 
netic energy density approaches a nonzero value for finite U. 
The incompatibility should indeed be there, as Phillips also 
realized in his next paragraph. Precisely, the quick rise of W\ 
from zero is the point which delivers the significance of the 
non-trivial finding. 

Beyond the quick rise, W\ appears to be consistent with 
a scaling 6 + rid which we proposed by considering the 
change in double occupancy rid- The overall behavior of W\ 
in this range of 5 can be reproduced by either our scaling 
W\ ~ 6 + rid or the perturbation result [6]. This is reason- 
ably understood by noting the fact that, as also mentioned in 
Ref. 0], the kinetic energy term in the perturbative expansion 
is closely related to the doublon density rid because the former 
creates doublon-holon pairs in the atomic limit. 

To explain the non-linear evolution of W\ at tiny doping 
followed by the scaling 5 + rid at higher doping, we need a 


new mechanism. It motivated us, in Ref. U], to propose an 
avalanchine screening for doublon-holon bound pairs. Once 
holes are doped into the Mott insulator whose gap is deter- 
mined by the doublon-holon binding, the mobile holes bring 
about the screening of the binding energy. The screened bind- 
ing energy survives as the pseudogap. What we proposed is 
that the screening should have a positive feedback accelerat- 
ing the weight transfer, since dissolved bound pairs further 
join in the screening process. This feedback will be particu- 
larly effective in 2D, because the screening is governed by the 
density of states at the Fermi level, which jumps to a nonzero 
value in 2D by doping with the formation of Fermi surface 
pocket. This self-accelerated process brings about the abrupt 
reduction of the doublon-holon binding and the non-linear in- 
crease of W\ upon doping. In Ref. HI, it was inferred to drive 
the system to the verge of a first-order transition or a marginal 
quantum criticality. 

In Comment J2l, Phillips claimed that we do not offer any 
explanation for the underlying mechanism of the quick reduc- 
tion of the Mott gap associated with W\ . On the contrary, as 
we repeated above, we already gave a picture in Ref. [1] which 
is qualitative but physically sound and, at least, reasonably fits 
to our numerical results. 

We reached, in Ref. [1], the conclusion that interplay be- 
tween doublon-holon bound pairs and mobile carriers is a key 
for understanding the Mott physics, as also addressed within a 
scenario of "2e-boson" recently hypothesized in Ref. |7[] and 
Comment 0]. To clarify the binding mechanism with quan- 
tum fluctuations and its variation as doping will give deeper 
understanding of the Mott physics beyond the prior perturba- 
tive arguments in the atomic limit. It would also be an intrigu- 
ing issue how our wide-ranging results including the hole- 
pockets, Fermi arcs, Lifshitz transitions and other nontrivial 
findings in Ref. [ 1 ] obtained by the state-of-the-art numerical 
calculations can be reproduced in the "2e-boson" calculations 
in Ref. 0]. 

Shiro Sakai, Yukitoshi Motome, and Masatoshi Imada 

Department of Applied Physics 

University of Tokyo 

Hongo, Tokyo 113-8656, Japan 


[1] S. Sakai, Y. Motome, and M. Imada, Phys. Rev. Lett. 102, 

056404 (2009). 

[2] P. Phillips, larXiv:0904.0454l 

[3] C. T. Chen, et al, Phys. Rev. Lett. 66, 104 (1991). 

[4] H. Eskes, M. B. J. Meinders, and G. A. Sawatzky, Phys. Rev. 

Lett. 67, 1035 (1991) 
[5] M. B. J. Meinders, H. Eskes, and G. A. Sawatzky, Phys. Rev. B 

48, 3916 (1993). 
[6] A. B. Harris and R. V. Lange, Phys. Rev. 157, 295 (1967). 
[7] P. Phillips, T.-P. Choy, and R. G. Leigh, Rep. Prog. Phys. 72, 

036501 (2009). 


